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Abstract 


The 6-year wind archives from the Goddard Institute for Space Stud- 
ies/Global Climate-Middle Atmosphere Model (GISS/GCMAM) were in- 
put to the GISS/Harvard/Irvine Chemical Transport Model (G/H/I CTM) 
to study the seasonal and interannual variability of the budgets and distri- 
butions of nitrous oxide (N2O) and trichlorofluoromethane (CCI3F), with 
the corresponding chemical loss frequencies recycled and boundary con- 
ditions kept unchanged from year to year. The effects of ozone feedback 
and quasi-biennial oscillation (QBO) were not included. However, the role 
of circulation variation in driving the lifetime variability is investigated. 
It was found that the global loss rates of these tracers are related to the 
extratropical planetary wave activity, which drives the tropical upward 
mass flux. For N2O, a semiannual signal in the loss rate variation is as- 
sociated with the interhemispheric asymmetry in the upper stratospheric 
wave activity. For CCI3F, the semiannual signal is weaker, associated 
with the comparatively uniform wave episodes in the lower stratosphere. 
The loss rates lag behind the wave activity by about 1-2 months. The 
interannual variation of the GCM generated winds drives the interannual 
variation of the annually averaged lifetime. The year-to-year variations 
of the annually averaged lifetimes can be about 3 % for N2O and 4 % for 
CCI3F. 


1 Introduction 

The relationship between the stratospheric wave activity and the meridional 
distribution of long-lived tracers was discussed by Holton (1986). The isopleths 
of a long-lived tracer are balanced by the transport due to diabatic circulation, 
horizontal eddy transport, and chemical loss in the stratosphere. The diabatic 
circulation steepens the isopleths, but both the horizontal eddy transport and 
the chemical loss tend to flatten the isopleths. Since both the diabatic circula- 
tion and the horizontal eddy transport are controlled by wave activity [Holton, 
1986; Holton et al., 1996], the effects of wave activity on the long-lived tracer 
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distributions are complicated. Greater wave activity can steepen the isopleths 
across subpolar latitudes by strengthening the diabatic circulation, but flatten 
the isopleths by enhancing the horizontal eddy transport. The effect of these 
opposing mechanisms was well illustrated by the different rates of descent of 
methane (CH4) and N 2 0 in the polar vortex [Strahan et al., 1996]. 

In the tropics, the upward mass flux across the tropopause is controlled by 
the extratropical wave activity [Haynes et ah, 1991; Holton et ah, 1996]. The 
annual mean of the mass flux across the tropical 100 mb surface is estimated to 
be 85x10 s kg/s with significant seasonal variation (114xl0 8 kg/s in December- 
January-February (DJF) and 56xl0 8 kg/s in June-July-August (JJA)) [Rosenlof 
and Holton, 1993]. This seasonal variation of the tropical upward mass flux in- 
duces an observable annual cycle of lower stratospheric temperatures in the 
tropics [Yulaeva et ah, 1994; Rosenlof, 1995 and references therein]. Interan- 
nual variation of the tropical upward mass flux was also derived from United 
Kingdom Meteorological Office (UKMO) assimilated temperatures and National 
Meteorology Center (NMC) temperatures [Rosenlof, 1995]. 

Both N 2 0 and CC1 3 F are long-lived tracers with tropospheric sources, and 
not effectively removed in the troposphere. They are transported into the strato- 
sphere through the tropical tropopause and dissociated in the stratosphere by 
chemical processes, providing constituents that trigger the catalytic destruc- 
tions of stratospheric ozone. The N 2 0 lifetime derived from observations is 
about 120 yrs [Prather et ah, 1994; Minschwaner et ah, 1998] and 124±49 yrs 
[Volk et ah, 1997], and that derived from previous 2-D models is in the range 
of 120-139 yrs [Ko and Jackman, 1994]. The CCI3F lifetime derived from ob- 
servations is 50 yrs [Prather et ah, 1994], 34-47 yrs [Ko and Jackman, 1994], 
52^ y rs [Cunnold et ah, 1997] and 41±12 yrs [Volk et ah, 1997], and that 
derived from the 2-D models is 40-60 yrs [Ko and Jackman, 1994]. 

In this paper, the seasonal and interannual variability in the N 2 0 and CCI3F 
lifetimes is examined with a 3-D model, which is required for investigation of 
the impacts of planetary wave on the tracers’ budgets since parameterization for 
the eddy transport are not necessary in 3-D models. The interannual variation 
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in ozone transport and quasi-biennial oscillation (QBO) are not included. This 
paper is the first step to understand the role of the variation in wave activity in 
driving the lifetime variations, without any feedback from ozone response and 
modulation by QBO. In section 2, the models used are described. Section 3 
discusses the simulation results. Discussion and conclusions are included in 
section 4. 


2 The Models 

The tracer experiments used the 21-layer Chemical Transport Model (CTM) 
[Prather, 1993], which advects chemical constituents with winds produced from 
the 23-layer GISS Global Climate/Middle Atmospheric Model (GISS/GCMAM) 
[Rind et ah, 1988a and 1988b]. The GISS/GCMAM has been analyzed for a 
range of climate change experiments [Rind et ah, 1990, 1992; Rind and Bal- 
achandran, 1995; Rind and Lonergan, 1995; Shindell et al., 1997], etc. The 
CTM adopts the 7.83°xl0° grid of the parent GCM for its horizontal resolution 
and uses vertical layers extending from the surface to 0.004 mb (about 85 km). 
The tropospheric layers, from 1 to 9, use sigma coordinates, while the strato- 
spheric layers, from 10 to 21, use fixed pressure levels beginning above 100 mb. 
The top three layers of the 23-layer GCM are combined into the top layer of the 
CTM. 

Boundary conditions for trace gas mixing ratios are applied in the lowest 
three layers (up to 1 km), which can be either fixed or varying in time. The 
chemistry of N2O and CCI3F is parameterized by their zonally and monthly 
averaged loss frequencies. These parameters are calculated from a photochemi- 
cal box model using observed climatic values for temperature, ozone, and other 
trace gases [Avallone and Prather, 1997; Prather and Jaffe, 1990]. Parameter- 
ized convection and the associated horizontal diffusion occur in the troposphere, 
but not in the stratosphere. The advecting scheme conserves first- and second- 
order moments of the tracer distribution [Prather, 1986]. The CTM has been 
used to study tropospheric distribution of CFCs, dilution of the ozone hole, the 
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space shuttle impacts, seasonal evolutions of N 2 0, O3, and C0 2 , and tracer- 
tracer correlations [Avallone and Prather, 1997; Hall and Prather, 1993, 1995; 
Prather et ah, 1987, 1990, 1993]. 

The GISS/GCMAM was used to produce 6-year wind archives and the as- 
sociated diagnostics, especially E-P flux convergence which is a representative 
for planetary wave activity. The model provides generally realistic interannual 
variability in the mean circulations and planetary wave activity [Rind et ah, 
1988b], The primary deficiencies of the model are somewhat reduced long- wave 
energy in the troposphere and lower stratosphere, too cold temperatures near 
the model top, and too warm temperatures in the Southern Hemisphere polar 
lower stratosphere. 

Initial tracer distributions were calculated by recycling the first year of the 
GCM winds until a near steady-state repeating annual cycle was reached. The 
subsequent 5 years of winds are then put into the CTM sequentially. These runs 
are referred to as SEQUENTIAL runs. In order to compare our results with 
the steady state results, we ran another set of experiments in which each year’s 
wind was recycled for about 10 years until a near steady-state was reached for 
that year’s wind archive. These are referred to as STEADYSTATE runs (i.e., 
w r e have 6 individual near steady-state distributions for the corresponding 6- 
year wind archives). Since we initialized the tracers with the first year steady 
state distributions, all the first year results in the SEQUENTIAL runs match 
those in the STEADYSTATE runs. For all the experiments, the chemical loss 
frequencies and the tropospheric convection processes (not very important in 
these experiments) were recycled from year to year. The fixed chemical loss fre- 
quencies are equivalent to fixing the overhead column ozone amount; therefore, 
feedbacks from ozone response are not allowed in this study. The boundary 
condition of N 2 0 was fixed at 300 ppbv, and that of CCI3F, at 250 pptv, for 
the lowest three levels of the CTM. 
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3 Results 


3.1 The Lifetimes 

The lifetime of a chemical tracer is defined as the total mass (in kg) divided 
by the globally integrated loss rate (in kg/yr). Since the boundary conditions 
were fixed, the total mass of both tracers did not vary significantly over the 
years. Hence, the lifetimes are inversely proportional to the loss rates in our 
simulation results. The difference in global loss rates between the SEQUEN- 
TIAL and STEADYSTATE runs is small (about 3-4% from the mean values, 
see figs. 5 and 6 ) for our 6 years of meteorology. The differences in tracer dis- 
tributions throughout the stratosphere are also small. Hereafter, our discussion 
will concentrate on the SEQUENTIAL runs. 

Table 1 shows the seasonal and interannual variability of the lifetimes of both 
tracers in the SEQUENTIAL runs. The largest interannual difference in annual- 
mean lifetimes is about 3% for N 2 0 and about 4% for CCI 3 F about their mean 
lifetimes. For both tracers, the seasonal variability is larger than the interannual 
variability. The interannual variability in the spring seasons (i.e. April (APR) 
or October (OCT)) is larger than that in the solstices (i.e. January (JAN) or 
July(JUL)). However, the interannual variability is reduced when the annual 
mean lifetime is taken. The model N 2 0 lifetime covers the range of 107-130 yrs, 
with the climatological mean of 121 yrs, consistent with those in Mahlman et al. 
(1986), Minschwaner et al. (1998), and Prather et al. (1994). The model CCI 3 F 
lifetime covers the range of 30-34 yrs, with the climatological mean of 31.3 yrs, 
lower than those in Cunnold et al. (1997), Ko and Jackman et al. (1994), and 
Volk et al. (1997). This discrepancy in the CCI 3 F lifetime may be due to the 
coarse vertical resolution (about 5 km in the stratosphere) applied by the CTM, 
which produces vertical mixing that is too rapid up to about 25 km [Avallone 
and Prather, 1997]. 
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3.2 The Distributions 


Since the N 2 0 loss is dominant in the upper stratosphere, the evolution of 
the zonal mean N 2 0 mixing ratio, averaged from 10 mb to 1 mb, is shown in 
figure la. In the tropics, there are primary peaks of mixing ratio, centered at 
about 10°S, during January to March, and secondary peaks of mixing ratio, 
centered at about 10°N, during June to September. These peaks of mixing 
ratio also correspond to the peaks of the loss rate shown in figure 5c. The 
location and time of the peaks agree with previous studies of the Cryogenic 
Limb Array Etalon Spectrometer (CLAES) data [Randel et ah, 1994; Lou et al., 
1997], in which primary peaks of N 2 0 mixing ratio were centered around 10°S 
during January to March in 1992 and 1993, while a secondary peak was centered 
around 10° N during August to October in 1992. However, it is hard to tell if the 
primary and secondary peaks had difference in magnitudes from these previous 
studies, although the model result of Lou et al. (1997) showed that the primary 
peaks had larger magnitude than the secondary peak at 10 mb. The annual 
cycle of the N 2 0 mixing ratio shown in figure la is also in general agreement 
with that of the CH 4 mixing ratio shown in Randel et al. (1998, their fig. 9). 
The primary peak of CH 4 mixing ratio at 2.2 mb was found during February to 
March, centered at about 10°S, and the secondary peak with about the same 
magnitude of the primary peak was found during August to September, centered 
at about 10°N. At 10 mb, the peak of CH 4 mixing ratio reached maximum at 
the equator from November to June [Randel et al., 1998; fig. 9 in the reference]; 
therefore, one can expect the secondary peak is smaller in magnitude if the CH 4 
mixing ratio is averaged from 10 mb to 1 mb. 

Figure lb shows the evolution of the zonal mean N 2 0 flux across 10 mb. In 
the tropics, primary peaks of upward flux occur in December- January, whereas 
secondary peaks of upward flux occur in June, with magnitudes smaller than 
those of the primary peaks. The evolution of the mixing ratio is consistent 
with that of the flux, with a time lag of 1-2 months, that is from mid-winter, 
when the flux is maximum, to late winter and early spring, when the mixing 
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ratio is maximum. In high latitudes, peaks of upward flux can be seen around 
70°N or 60°S, where polar descent (downward motion) should be expected. The 
zonal mean tracer flux includes the mean and eddy fluxes, and resembles the 
pattern of zonal mean circulation, instead of the Transformed Eulerian Mean 
(TEM or the residual) circulation, because the eddy flux is about 2-3 order 
of magnitude smaller. Therefore, the pattern of the zonal mean flux shown 
in figure lb (as well as in fig. 2 b for CCI 3 F) is consistent with the multi-cell 
structure of the zonal mean circulation, and the upward flux seen around 70°N 
or 60° S is an extension of the upward branch of the polar and the Ferrel cells 
from the troposphere. However, the polar upward flux occurring during polar 
nights less affects the lifetime and its variation because the loss is dominant in 
the tropics to the subtropics. 

Similar results can be obtained for CCI 3 F if we consider a region in the lower 
stratosphere, because of the predominant CCI 3 F loss in this region. Figure 2a 
shows the zonal mean CCI 3 F mixing ratio, averaged from 100 mb to 10 mb, 
and figure 2b, the CCI 3 F upward flux across 100 mb. The seasonal evolution 
of subtropical to mid-latitude CCI 3 F mixing ratios (see the contour lines of 
105-150 pptv in fig. 2a) agree with that of the CH 4 mixing ratio at 68 mb in 
Randel et al. (1998), who found that the contour lines shift northward during 
July and southward during December to January as in our case of CCI 3 F mixing 
ratios. However, the CCI 3 F mixing ratio in the tropics shows less variability 
than those of N 2 0 and CH 4 at the same levels. This may be due to the fact 
that CCI 3 F has larger photochemical loss in the tropics and behaves less as an 
passive tracer in this region. The seasonal variability of CCI 3 F mixing ratio 
recovers in higher latitudes where the photochemical loss is smaller. 

3.3 Planetary Wave Activity and Upward Tracer Fluxes 

The extratropical ’’suction pump” is believed responsible for driving the tropical 
upward mass flux which brings in tracers from the troposphere [Haynes et al., 
1991; Holton et al., 1996]. The theory suggests that the tropical upward mass 
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flux is related to the extratropical wave activity, including both planetary and 
subgrid scale waves. Here we will consider only the impacts of planetary wave 
activity. Figure 3 shows the model evolution of planetary wave activity, repre- 
sented by the E-P flux convergence averaged from 10 mb to 1 mb in figure 3a, and 
from 100 mb to 5 mb in figure 3b. In order to justify that the GISS/GCMAM 
has proper variability in planetary wave activity, the E-P flux divergence de- 
rived from NMC data from 1981-1995 is shown in figure 4 [Randel, 1992] (note 
that fig. 4 plots the E-P flux divergence so the shaded negative regions should 
be compared with the red positive regions in fig. 3). In the upper stratosphere 
(10-1 mb, fig. 3a and top panel of fig. 4), the model E-P flux convergence in 
middle latitudes reaches 2-3 m/s/day in northern hemisphere, and 1-2 m/s /day 
( 2-3 m/s/day in the 4-th year) in southern hemisphere. This is consistent with 
the NMC E-P flux data, about 2.5 m/s/day in most of the years (5 m/s/day 
in some of the years) in northern hemisphere, and about 0-2.5 m/s/day (over 
2.5 m/s/day in some of the years) in southern hemisphere. In the northern 
polar region, maxima over 5 m/s/day are seen in some of the years in both the 
model and NMC data. In the lower stratosphere (100-5 mb, fig. 3b and bottom 
panel of fig. 4), the model E-P flux convergence generally agrees with the NMC 
E-P flux convergence (2-3 m/s/day in northern hemisphere and 1-2 m/s/day 
in southern hemisphere). In both upper and lower stratosphere, the model can 
provide reasonable variability in planetary wave activity. 

In the upper stratosphere (fig. 3a), strong wave episodes occur during au- 
tumn to early spring, from 8 ° to 40° in latitude, and drive larger tropical upward 
tracer flux. The interhemispheric asymmetry in the winter wave activity is as- 
sociated with the seasonal asymmetry in the tropical upward flux, as well as 
the seasonal asymmetry in the N 2 0 loss. In the lower stratosphere (fig. 3b), 
wave episodes cover a range in time and latitude (from 8 ° to 72°) that is wider 
than those in the upper stratosphere. The interhemispheric asymmetry in wave 
activity is less evident in the lower stratosphere. The distinction between pri- 
mary and secondary peaks of the tropical CCI 3 F flux and mixing ratio is also 
less evident. 
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The E-P flux convergence shown in Figure 3 is consistent with the tracer 
flux patterns shown in Figures lb and 2b. In the upper stratosphere, peaks of 
downward tracer flux across 10 mb (fig. lb) occur when the tropical upward 
tracer flux is at its maxima. This is because the E-P flux convergence from 
20° to 40° (fig. 3a) provides poleward transport while the E-P flux divergence 
from 50° to 70° provides equatorward transport. The resultant horizontal tracer 
flux convergence produces downward tracer flux from 40° to 50° due to mass 
continuity. The downward tracer flux continues to lower altitude as shown in 
Figure 2b. In the 100 mb level (fig. 2b), the upward tracer flux appearing in 
the southern polar region is stronger than that in the northern polar region 
and may shift the spring time total ozone maxima toward lower latitudes in the 
southern hemisphere. This is consistent with the previous observational and 
model evolution of the total ozone, in which maxima of total ozone are seen 
during southern spring in the southern subpolar region (about 400 DU around 
50-60°S in September to October), whereas the maxima of total ozone during 
northern spring occur within the northern polar region (about 440 DU around 
80-90°N in March to April) [Bowman and Krueger, 1985; Ushimaru and Tanaka, 
1994]. 

3.4 The Time Series 

Time series of several averaged variables, representing different physical pro- 
cesses, are shown in figures 5 and 6 for N2O and CCI3F, respectively. Different 
tracers have different regions of dominant loss, and so only data in appropriate 
regions are used to produce the corresponding time series. Figure 5a is the 
averaged E-P flux convergence from 8° to 40° in latitude and from 10 mb to 
1 mb in altitude, a region where the planetary wave activity is responsible for 
the tropical upward flux across 10 mb. The dotted and dashed lines represent 
the northern and southern wave activity, respectively. The tropical upward flux 
is driven by the combination of wave activity in both hemispheres; therefore, 
the sum is also plotted as the solid line. The spikes occurring during November 
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to March are mainly from the northern wave episodes, whereas those occurring 
during May to October are mainly from the southern wave episodes. Figure 5b 
is the N 2 0 flux averaged from 24°S to 24°N, a latitude range in which the N2O 
flux is upward in the tropics, at 10 mb (fig. lb) . The time series of the flux is 
correlated with that of the wave activity. 

Figures 5c and d show the variations of the global N2O loss rate and mixing 
ratio, respectively. The mixing ratio is averaged from 30°S to 30°N in latitude 
and from 10 mb to 1 mb in altitude, a region in which the averaged N 2 0 mixing 
ratio reflects the tropical concentration variation (see fig. la). The solid lines arc 
from the STEADYSTATE runs and the dotted lines are from the SEQUENTIAL 
runs. Primary peaks of mixing ratio occur during January to March, whereas 
secondary peaks occur during August to October, with magnitudes smaller than 
those of the primary peaks. This semiannual signal in the loss rate and the 
mixing ratio resembles the semiannual signal in the tropical upward flux across 
10 mb (fig. 5b), which is driven by the asymmetric wave forcing from both 
hemispheres (fig. 5a). This result is consistent with Rosenlof (1995) in which a 
semiannual signal in the tropical upward mass flux above 10 mb was detected 
(fig. 15 in the reference). 

The CCI3F loss is dominant in the lower stratosphere, over a range of lati- 
tude wider than that found for the N 2 0 loss. Figure 6a is the averaged E-P flux 
convergence from 8° to 72° in latitude and from 100 mb to 5 mb in altitude, 
a region in which the averaged E-P flux convergence represents about 80% of 
the forcing of the tropical upward mass flux [Rosenlof, 1995]. The symbols in 
figure 6a are similar to those in figure 5a. Peaks of wave episodes are compar- 
atively uniform in time. However, primary peaks during January to April are 
still evident. Figure 6b shows the time series of the upward CCI3F flux across 
the tropical tropopause, represented by the flux averaged from 30° S to 30° N at 
100 mb, a region where the tropical flux is mainly upward. The annual cycle is 
evident with peaks occurring during January to March. The time series of the 
CCI3F loss rate and mixing ratio are shown in Figure 6c and d, respectively. 
The mixing ratio is averaged from 50°S to 50°N in latitude and from 100 mb 
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to 10 mb in altitude, a region with wider latitude range because the dominant 
loss of CCI3F covers a wider range in latitude. The semiannual signal is less 
evident than that in the N2O plots in figure 5, associated with the compara- 
tively uniform wave episodes in the lower stratosphere. This is consistent with 
Rosenlof (1995), who found that the semiannual signal has a larger amplitude 
above 10 mb. 

An estimate of the air mass flux across 100 mb can be obtained if the tracer 
fluxes and the corresponding tracer mixing ratios right below 100 mb are known. 
Table 2 shows the estimate of the tropical upward air mass flux across 100 mb 
from the CTM as well as the estimates from the previous literature. Our air mass 
flux is smaller than the estimate by Rosenlof and Holton (1993) during northern 
winter. However, during southern winter, our air mass flux is quite consistent 
with Rosenlof and Holton (1993). This is reasonable because GISS/GCMAM 
simulates weaker long-wave activity and the coarse resolution cannot resolve 
smaller-scale processes which may govern the troposphere-stratosphere exchange 
[Rind et al., 1988a]. Also shown in table 2 is the estimate by Rosenlof (1995) 
at 70 mb using the UKMO assimilated temperatures. Our CTM is consistent 
with Rosenlof (1995), recalling that our data were sampled at a low T er altitude 
(100 mb) and thus should give slightly larger values than those at 70 mb. For the 
interannual variability in the tropical upward tracer fluxes (fig. 5b and 6b), our 
CTM also gives reasonable variability w T hich is similar to that showm in Rosenlof 
(1995) (their figure 11 shows a 2-year time series of the tropical upward mass 
flux at 70 mb). 

3.5 The Time-Lag 

From the correlations of the time series shown in figures 5 and 6, one can 
obtain the time lag between different physical processes. Figure 7 shows these 
cross correlations between the adjacent time series, which are obtained from the 
following formula: 
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where m is the time lag from -12 to 12 months, overbar variables represent the 
mean (the average of 72 monthly data except for wave activity). For the wave 
activity, the mean is chosen to be the offset, which is zero in figure 5a but 1.0 
in figure 6a, because the wave activity appears as spikes instead of a Gaus- 
sian distribution around the mean. This choice of the means for wave activity 
increases the correlation but does not change the time-lag between different pro- 
cesses. Since there are clear annual cycles in the time series, peaks of correlation 
around ±12 months and anticorrelation around ±6 months are expected. For 
N 2 0 (fig. 7a), the tropical upward flux at 10 mb lags the upper stratospheric 
wave activity by 1 month (dotted line); the loss lags the tropical upward flux 
by less than 1 month (dashed line); and the tropical upper stratospheric mixing 
ratio lags the loss by about 1 month (dash-dot line). The direct result is that the 
mixing ratio lags the wave activity by about 2 months (solid line) . This result 
is consistent with Randel et al. (1998) who found a time-lag about 1 month 
between the CH4 maxima at 2.2 mb and the tropical upward motion at the same 
level. Since the horizontal eddy transport counteracts the transport caused by 
the mean vertical wind in the tropics, the years of the maximum wave activity 
and the maximum tropical upward flux does not coincide for N 2 0 (the second 
year for the wave activity but the third year for the tropical upward flux). 

The story is different for CCI3F (fig. 7b) because its predominant loss is in the 
lower-middle stratosphere. The correlation between the tropical upward flux and 
the lower stratospheric wave activity is small (about 0.3) and covers a broader 
time range (dotted line), which is associated with the comparatively uniform 
wave episodes in the lower stratosphere. The loss responds quite instantaneously 
to the tropical upward flux (time lag is less than 1 month, dashed line); and the 
lower stratospheric mixing ratio lags the loss by less than 1 month (dash-dot 
line). The direct correlation between the mixing ratio and the wave activity is, 
as expected, widespread, and it peaks at about 1 month time lag. 
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In our work, only planetary wave activity has been considered for the cor- 
relation calculations. GISS/GCMAM has parameterized gravity wave drag 
[Rind et al., 1988a] which can also drive the tropical upward mass flux. In- 
clusion of the gravity wave drag into the total wave activity for the correlation 
calculations (not shown) changes the peaks of correlation by about 10-20% but 
does not change the time lags discussed above. 

4 Discussion and Conclusions 

According to the downward control mechanism, the extratropical wave activity 
drives the tropical upward mass flux, which brings long-lived tracers with tro- 
pospheric sources into the stratosphere. The stratospheric tracer distribution is 
determined by a balance between transport and photochemical loss. For N2O, 
the loss is dominant in the upper stratosphere, and its upper stratospheric mix- 
ing ratio lags the upper stratospheric wave activity by about 2 months. Strong 
wave episodes in winter induce peaks of mixing ratio in late winter to early 
spring. The interhemispheric asymmetry in wave activity is associated with 
a clear semiannual signal in the seasonal variation of the N2O mixing ratio. 
For CCI3F, the loss is dominant in the lower-middle stratosphere, where wave 
episodes are comparatively uniform in time although northern winter peaks are 
still evident. The resultant seasonal variation of the CCI3F mixing ratio has an 
evident annual cycle but a vague semiannual signal. 

QBO, which was not simulated in our models, modulates the tropical upward 
mass flux in the stratosphere and contributes the majority of the interannual 
variance in stratospheric tracer concentrations (e.g., upper stratospheric CH 4 
investigated in Randel et al. (1998)). Our model studies here cannot show the 
impacts of QBO on the interannual variation of the tracer budgets. However, 
the role of the model generated variation in planetary wave activity can then 
be studied without the modulation of QBO. Inclusion of QBO and the relative 
contributions of different natural forcings (e.g., QBO, ENSO, and aerosol effects) 
in driving the interannual variation in the tropical upward mass flux need further 
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studies by both model and observational data analyses. 

In the literature, people working on CTMs usually recycle a 1-year wind 
archive until the models reach steady states. In this study, we see that the dif- 
ference between these steady-state runs and runs using subsequent wind archives 
with interannual variation should be small. The wind archives generated by the 
GISS/GCMAM produce clear seasonal and interannual variations in both the 
N 2 0 and CCI3F budgets (e.g., their stratospheric concentrations and loss rates). 
This illustrates how climate change (here we have the GCM generated wind vari- 
ability) can affect long-lived tracer budgets through dynamical processes. The 
variation of the model tropical upward flux across 100 mb is consistent with the 
air mass flux variation obtained by the radiatively determined residual circu- 
lation [Rosenlof, 1995]. The magnitude of the CTM air mass flux in northern 
winter is smaller than the result in Rosenlof and Holton (1993), but that in 
southern winter is consistent with them. 

Comparing the tropical upward tracer fluxes across 100 mb (not shown for 
N 2 0 and fig. 6b for CCI3F) with the total loss rates (fig. 5c and 6c for N 2 0 
and CCI3F, respectively), we find that about one-tenth of N 2 0 and one-third of 
CCI3F entering the stratosphere from the tropical tropopause is dissociated in 
the stratosphere (the ratio is about the ratio of the lifetimes of the two tracers if 
the tropospheric mixing ratios are fixed as they are in our case). The remaining 
amount will be transported back to the troposphere in middle to high latitudes. 
The variation of the ventilation through tropical tropopause may affect the 
concentrations of other tracers in the upper troposphere, implying variations of 
the upper tropospheric chemistry in the tropics. Moreover, the variation of the 
mid-latitude downward flux across 100 mb implies variation of the extratropical 
stratosphere-troposphere exchange. 

The annual mean lifetimes show the largest year-to-year variability of about 
3-4% about the mean values of 121 yrs for N 2 0 and 31 yrs for CCI3F, smaller 
than the seasonal variability which is over 10% about the mean values. Upper 
Atmospheric Research Satellite (UARS) may provide data necessary to verify 
the model variations in the stratospheric mixing ratios and the tracer lifetimes 
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[Minschwaner et al., 1998]. 

The theoretical framework by Plumb and Ko (1992) suggests a relation be- 
tween the ratio of two tracers’ lifetimes and the slope of the tracer-tracer cor- 
relation plot in the lower stratosphere, 


n ~ 

72 02 dtp l 

where r is the tracer lifetime, 0 is the tracer mixing ratio, and represents 
the slope of the tracer-tracer correlation plot. Scatter of data points is expected 
if both tracers undergo different chemistry in the stratosphere [Avallone and 
Prather, 1997; Hall and Prather, 1995]. Using equation (2), one can relate the 
interannual variability of tracer-tracer correlation to the interannual variability 
of tracer lifetimes, that is, 


#2 #2 ^ An Ar 2 

# i n 72 


( 3 ) 


where we have assumed the tropospheric mixing ratios do not have an interan- 
nual variation as they do in our simulations. From table 1 and equation (3), 
one can imply that the interannual variability of the CCI 3 F and N 2 O correla- 
tion plot is about 1-4%. Therefore, one cannot expect to detect it by looking 
at the model or observational tracer-tracer correlations because the signal is 
far smaller than the scatter of the plot [Avallone and Prather, 1997]. One can 
also infer that the variability of lifetimes discussed above has a very small effect 
on the estimates of several environmental indices which depend on the ratio of 
lifetimes of the tracers being concerned (e.g., Chlorine Loading Potential, etc). 
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Figure 1. Evolution of the N2O distribution and flux, (a) The zonal mean 
N 2 0 mixing ratio (in ppbv) averaged from 10 mb to 1 mb; (b) the zonal mean 
N 2 0 flux (in 10 -12 kg/m 2 /s) across 10 mb. The x-axis is for simulation year 
and the y-axis is for latitude. 

Figure 2. Evolution of the CCI3F distribution and flux, (a) The zonal 
mean CCI3F mixing ratio (in pptv) averaged from 100 mb to 10 mb; (b) the 
zonal mean CCI3F flux (in 10“ 13 kg/m 2 /s) across 100 mb. The x-axis is for 
simulation year and the y-axis is for latitude. 

Figure 3. Evolution of the model stratospheric planetary wave activity, 
represented by the E-P flux convergence (in m/s/day), (a) averaged from 10 mb 
to 1 mb; (b) averaged from 100 mb to 5 mb. The x-axis is for simulation year 
and the y-axis is for latitude. 

Figure 4. Evolution of the NMC stratospheric planetary wave activity, 
represented by the E-P flux divergence (in m/s/day), averaged from 10 mb to 
1 mb (top panel, contour interval 2.5), and from 100 mb to 5 mb (bottom panel, 
contour interval 0.5) [Randel, 1992]. 

Figure 5. Time series of several averaged variables for N 2 0. (a) The E-P 
flux convergence (in m/s/day) averaged from 8° to 40° in latitude and from 
10 mb to 1 mb in altitude, with the dotted line representing the northern hemi- 
sphere, the dashed line the southern hemisphere, and the solid line the sum 
of values from both hemispheres; (b) the tropical upward N 2 0 flux (in 10“ 12 
kg/m 2 /s) across 10 mb, averaged from 24° S to 24°N in latitude; (c) the global 
N 2 0 loss rate (in 10“ 12 kg/m 2 /s); (d) the zonal mean N 2 0 mixing ratio (in 
ppbv) averaged from 30°S to 30°N in latitude and from 10 mb to 1 mb in 
altitude. 
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Figure 6. Time series of several averaged variables for CCI3F. (a) The 
E-P flux convergence (in m/s/day) averaged from 8° to 72° in latitude and 
from 100 mb to 5 mb in altitude, with the dotted line representing the northern 
hemisphere, the dashed line the southern hemisphere, and the solid line the sum 
of values from both hemispheres; (b) the tropical upward CCI3F flux (in 10 -14 
kg/m 2 /s) across 100 mb, averaged from 30°S to 30°N; (c) the global CCI3F 
loss rate (in 10 -14 kg/m 2 /s); (d) the zonal mean CCI3F mixing ratio (in pptv) 
averaged from 55°S to 55°N in latitude and from 100 mb to 10 mb in altitude. 

Figure 7. Time-lag correlations for (a) the adjacent N2O time series from 
Figure 4; (b) the adjacent CCI3F time series from Figure 5. Dotted lines repre- 
sent correlations between the extratropical wave activity and the upward flux, 
dashed lines, between the upward flux and the global loss rate, dash-dot lines, 
between the global loss rate and the averaged mixing ratio, and solid lines, 
between the wave activity and the averaged mixing ratio. 
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Table 1. The seasonal and interannual variability in the N 2 0 and CCI 3 F 
lifetimes (in yrs). 



JAN 

APR 

JUL 

OCT 

Annual mean 

n 2 o 

1 . 

109.3 

118.4 

130.2 

129.8 

122.2 

2. 

109.8 

120.5 

125.5 

124.8 

120.2 

3. 

109.6 

118.6 

129.7 

129.8 

121.1 

4. 

108.5 

125.1 

129.7 

128.4 

123.2 

5. 

107.2 

115.7 

129.6 

125.4 

119.9 

6. 

110.7 

119.0 

126.3 

123.6 

120.5 

Average 

109.2 

119.6 

128.5 

127.0 

121.2 

Max. Changes (%) 

3.2 

7.9 

3.7 

3.9 

2.7 

CCI 3 F 

1 . 

30.3 

30.1 

34.1 

31.2 

31.5 

2. 

29.8 

29.7 

33.2 

30.0 

30.6 

3 . 

29.9 

30.5 

34.0 

31.4 

31.3 

4. 

30.2 

31.3 

34.3 

31.1 

31.8 

5 . 

29.7 

30.3 

34.2 

30.7 

31.5 

6. 

31.1 

31.1 

33.6 

29.0 

31.1 

Average 

30.2 

30.5 

33.9 

30.6 

31.3 

Max. Changes (%) 

4.6 

5.2 

3.2 

7.8 

3.8 


Note: In last column are the annual mean lifetimes, and in the last row of 
each section are the largest peak-to-peak interannual differences in the corre- 
sponding columns (in percentage about the corresponding column mean values) . 
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Table 2. Comparison of the tropical upward mass flux (in 10 8 kg/s). 


Data Source 

DJF 

JJA 

Annual mean 

Rosenlof and Holton (1993) 

114 

56 

85 

Rosenlof (1995) (at 70 mb) 

66.4 

43.2 

54.8 

Our CTM 

88 

50 

69 


Note: The data of Rosenlof (1995) are presented at 70 mb and for January 
and July only, instead of for the DJF and JJA seasons. The last column is 
simply the average of the DJF and JJA column data. 
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